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In this work, we find exact gravastar solutions in the context of noncommutative geometry, and 
explore their physical properties and characteristics. The energy density of these geometries is 
a smeared and particle-like gravitational source, where the mass is diffused throughout a region 
of linear dimension ^/a due to the intrinsic uncertainty encoded in the coordinate commutator. 
These solutions are then matched to an exterior Schwarzschild spacetime. We further explore the 
dynamical stability of the transition layer of these gravastars, for the specific case of /3 = KP/a < 1.9, 
where M is the black hole mass, to linearized spherically symmetric radial perturbations about 
static equilibrium solutions. It is found that large stability regions exist and, in particular, located 
sufficiently close to where the event horizon is expected to form. 
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I. INTRODUCTION 



Recently, an alternative picture for the final state of 
gravitational collapse has emerged [1|]. The latter, de- 
noted as a gravastar ((^razitational uacuum star), con- 
sists of an interior compact object matched to an exte- 
rior Schwarzschild vacuum spacetime, at or near w^here 
the event horizon is expected to form. Therefore, these 
alternative models do not possess a singularity at the ori- 
gin and have no event horizon, as its rigid surface is lo- 
cated at a radius slightly greater than the Schwarzschild 
radius. More specificallvjthe gravastar picture, proposed 
by Mazur and Mottola fl| , has an effective phase transi- 
tion at / near where the event horizon is expected to form, 
and the interior is replaced by a de Sitter condensate. 
This new emerging picture consisting of a compact ob- 
ject resembling ordinary spacetime, in which the vacuum 
energy is much larger than the cosmological vacuum en- 
ergy, is also denoted as a "dark energy star" 0. In fact, 
a wide variety of gravastar models have been considered 
in the literature [3, |3| and their observational signatures 
have also been explored [Ij. In this work, we consider a 
further extension of the gravastar picture in the context 
of noncommutative geometry. The dynamical stability 
of the transition layer of these gravastars to linearized 
spherically symmetric radial perturbations about static 
equilibrium solutions is also explored. The analysis of 
thin shells [Q] and the respective linearized stability anal- 
ysis of thin shells has been recently extensively considered 
in the literature, and we refer the reader to Refs. 0, ll| 



for details. 

In the context of noncommutative geometry, an in- 
teresting development of string/M-theory has been the 
necessity for spacetime quantization, where the space- 
time coordinates become noncommuting operators on a 
ZJ-brane The noncommutativity of spacetime is en- 
coded in the commutator [x^, x"] — i 9^'^ , where 6"^" is an 
antisymmetric matrix which determines the fundamen- 
tal discretization of spacetime. It has also been shown 
that noncommutativity eliminates point-like structures 
in favor of smeared objects in flat spacetime [loj . Thus, 
one may consider the possibility that noncommutativity 
could cure the divergences that appear in general relativ- 
ity. The effect of the smearing is mathematically imple- 
mented with a substitution of the Dirac-delta function 
by a Gaussian distribution of minimal length ^/a. In 
particular, the energy density of a static and spherically 
symmetric, smeared and particle-like gravitational source 
has been considered in the following form [ll| 



M ( \ 



(1) 



where the mass M is diffused throughout a region of lin- 
ear dimension ^Ja due to the intrinsic uncertainty en- 
coded in the coordinate commutator. 

The Schwarzschild metric is modified when a non- 
commutative spacetime is taken into account [ill . IT^ . 
The solution obtained is described by the following space- 
time metric 



ds^ = -fir) dt^ 



dr^ 
W) 



+ (de^ + sin^ d(l)^) , (2) 
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with f{r) = 1 — 2m{r)/r, where the mass function is 
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defined as 



following relationships 



m(r) 



2M ^ 



(3) 



and 
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is the lower incomplete gamma function [ 1 1| . The classi- 
cal Schwarzschild mass is recovered in the limit r/^/a — > 
oo. It was shown that the coordinate noncommutativity 
cures the usual problems encountered in the description 
of the terminal phase of black hole evaporation. More 
specifically, it was found that the evaporation end-point 
is a zero temperature extremal black hole and there ex- 
ist a finite maximum temperature that a black hole can 
reach before cooling down to absolute zero. The existence 
of a regular de Sitter at the origin's neighborhood was 
also shown, implying the absence of a curvature singular- 
ity at the origin. Recently, further research on noncom- 
mutative black holes has been undertaken, with new solu- 
tions found providing smeared source terms for charged 
and higher dimensional cases p^ . Furthermore, exact 
solutions of semi-classical wormholes [3] in the context 
of noncommutative geometry were found and their 
physical properties and characteristics were analyzed. 

This paper is outlined in the following manner. In 
Section ini we present the generic structure equations of 
gravastars, and specify the mass function in the context 
of noncommutative geometry. In Section IIIIl the lin- 
earized stability analysis procedure is briefly outlined, 
and the stability regions of the transition layer of gravas- 
tars are determined. Finally in Section [IVI we conclude. 
We adopt the convention G = c = 1 throughout this 
work. 



II. STRUCTURE EQUATIONS OF 
GRAVASTARS IN NONCOMMUTATIVE 
GEOMETRY 



Consider the interior spacetime, without a loss of gen- 
erality, given by the following metric, in curvature coor- 
dinates 



1 — 2m(r)/r 



+ dn"^ , (5) 



where dfi^ = {dO"^ + sin^ 8 d(jp); $(r) and m(r) are arbi- 
trary functions of the radial coordinate, r. The function 
m(r) is the quasi-local mass, and is denoted as the mass 
function. 

The Einstein field equation, G^^ — SttT^^ provides the 



m 

P'r 



r{r — 2m) ' 
{p + Pr)im 



r{r — 2m) 



(6) 
(7) 

+ -(Pt-Pr), (8) 



where the prime denotes a derivative with respect to 
the radial coordinate. p{r) is the energy density, Priv) 
is the radial pressure, and pt{r) is the tangential pres- 
sure. Equation Q corresponds to the anisotropic pres- 
sure Tolman-Oppenheimer-Volkoff (TOY) equation. The 
factor <I>'(r) may be considered the "gravity profile" as 
it is related to the locally measured acceleration due 
to gravity, through the following relationship A = 
•\/l — 2m{r) / r ^' [r) . The convention used is that $'(r) 
is positive for an inwardly gravitational attraction, and 
negative for an outward gravitational repulsion. 

Using the equation of state, pr = — p, and taking into 
account the field equations ([6]) and ([7]), we have the fol- 
lowing relationship 



$'(r) 



r (r — 2m) 
which provides the solution given by 



$(r) = i In 
^ ' 2 



1 - 



2m(r) 



(9) 



(10) 



One now has at hand three equations, namely, the field 
Eqs. ©-([SI), with four unknown functions of r, i.e., p{r), 
Pr{r), ptir), and m{r). We shall consider the approach 
by choosing a specific choice for a physically reasonable 
mass function m(r), thus closing the system. 

In this context, we are interested in the noncommuta- 
tive geometry inspired mass function given by Eq. 
The latter is reorganized into the following form 



m(r) 



2M 
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where /3 is defined as /3 = AP/a. 

Note that three cases need to be analyzed [Tll |: 

a) If /3 < 1.9, no roots are present; 

b) If /3 > 1.9, we have two roots, r_ and r+, with r+ > 

c) If /3 = 1.9, we have r+ = r_, which may be inter- 

preted as an extreme situation, such as the extreme 
Reissner-Nordstrom metric. 

The function /(r) = (1 — 2m{r)/r) is depicted in Fig. 
[l]for these three cases for the following values /3 = 1.5, 
/3 = 1.9 and /3 = 3.5, respectively. Note that all the roots 
lie within the Schwarzschild radius = 2M, where M is 
the total mass of the system. 
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and ([12]), respectively, the four- velocity of the junction 
surface x^{t, 9, 0) = (^(t), a(T), 0, </>) is given by 



'1 - 



1 



-,a,0,0 ,(14) 



where the (±) superscripts correspond to the exterior and 
interior spacetimes, respectively, so that m± are defined 
as rn_ = m(a) and m+ = Af, respectively. 

The unit normal 4— vector, n^, to E is defined as 



= ± 



Q/3 



df df 



-1/2 



dl_ 

dxf' 



(15) 



FIG. 1: The function /(r) = (1 - 2m{r)/r) is depicted for 
the three cases with the following values /3 = 1.5, j3 = 1.9 and 
P = 3.5, respectively. 



III. LINEARIZED STABILITY OF 
GRAVASTARS IN NONCOMMUTATIVE 
GEOMETRY 



with rt^ = +1 and n^e'^.-^ = 0. The Israel formalism 
requires that the normals point from the interior space- 
time to the exterior spacetime [T6| . Thus, for the in- 
terior and exterior spacetimes given by the metrics ([S]) 
and p2|) . respectively, the normals may be determined 
from Eq. or from the contractions U^n^ — and 



-|-1, and are provided by 



A. Junction interface and surface stresses 

We shall model specific gravastar geometries by match- 
ing an interior gravastar geometry, given by Eq. ([5l) . 
where the metric functions are given by Eqs. (jlOp and 
([TT|). with an exterior Schwarzschild solution 



r I 



2M\ 



dr' 



[de"^ + SIT? e dcj)^) , 



(12) 



at a junction interface S, situated outside the event hori- 
zon, a > rb = 2M. We emphasize that the larger root 
r_|_ lies inside the Schwarzschild event horizon. Thus, in 
this work we are only interested in the case of /3 < 1.9, 
which corresponds to the absence of event horizons for 
the inner solution. 

Consider the junction surface E as a timelike hyper- 
surface defined by the parametric equation of the form 
fix^{^^)) =0. C = i'Tj (f) are the intrinsic coordinates 
on E, where r is the proper time on the hypersurface. 
The three basis vectors tangent to E are given by e^j) = 



9/9^*, with the following components 



dxt'/dC 



The induced metric on the junction surface is then pro- 
vided by the scalar product gij = e(i) • e^) = gf^^e'^.^e'^^y 
Thus, the intrinsic metric to S is given by 



-dr' 



{d0^ 



(13) 



Note that the junction surface, r = a, is situated outside 
the event horizon, i.e., a > ri,, to avoid a black hole 
solution, and we are only interested in the case of /3 < 1.9, 
of Eq. pT|) . as emphasized above. 

For the specific cases considered in this work, namely, 
the interior and exterior spacetimes given by Eqs. ([S]) 



-,0,0 



(16) 



respectively, with m± defined as to_ 
M, as before. 



TO (a) and toj 



The extrinsic curvature is defined as Kij = n^.^i,e^^-^e^^^-y 



Differentiating n^e^^ 



with respect to we have 



ture is finally given by 



I, ^r^r, so that the extrinsic curva- 



± dx^dx^\ 



(17) 



Note that, in general, Kij is not continuous across E, so 
that for notational convenience, the discontinuity in the 



extrinsic curvature is defined 



Taking into account the interior spacetime metric ([S]) 
and the Schwarzschild solution (jl2p . the non-trivial com- 
ponents of the extrinsic curvature are given by 



M 



2m(a) 



(18) 



(19) 



and 







1 - 



2M 



2m{a) 



(20) 
(21) 
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respectively. 

The Einstein equations may be written in the following 
form, 



1 



(22) 



denoted as the Lanczos equations, where 5^ is the sur- 
face stress-energy tensor on S. Considerable simplifi- 
cations occur due to spherical symmetry, namely k*^- = 
diag K^g, K%). The surface stress-energy tensor may 
be written in terms of the surface energy density, a, and 
the surface pressure, V, asSj = diag(— cr, P, P). Thus, 
the Lanczos equation, Eq. (I22L then provide us with the 
following expressions for the surface stresses 



V = 




1 - 



2m 



n2 



(23) 



1 _ ™ _ to' + a. 



1 - 



We also use the conservation identity given by S[ 



0) 



where \X\^ denotes the discontinuity 

across the surface interface, i.e., \X\^ — X+js — X^js. 
The momentum fiux term in the right hand side corre- 
sponds to the net discontinuity in the momentum fiux 
— T/ji/ which impinges on the shell. The conserva- 
tion identity is a statement that all energy and momen- 
tum that plunges into the thin shell, gets caught by the 
latter and converts into conserved energy and momentum 
of the surface stresses of the junction. 

Note that S'^jj = — [tr + 2d((T + 7')/a], so that the con- 
servation identity provides us with 



a' = {o + V). 

a 



(25) 



This relationship will be used in the linearized stability 
analysis considered below. 



B. Linearized stability analysis 

Using the surface mass of the thin shell nis = Ana'^a, 
Eq. (j25p can be rearranged to provide the following re- 
lationship 



(26) 



with the parameter i] defined as i] — V' /cr' , and T given 
by 



An 

T = — (a + 7') 
a 



Equation (j26p will play a fundamental role in deter- 
mining the stability regions of the respective solutions. 



Note that 77 is used as a parametrization of the stable 
equilibrium, so that there is no need to specify a surface 
equation of state. The parameter is normally inter- 
preted as the speed of sound, so that one would expect 
that < r] < 1, based on the requirement that the speed 
of sound should not exceed the speed of light. We re- 
fer the reader to Refs. [8] for further discussions on the 
respective physical interpretation of 77 lying outside the 
range < 77 < 1. 

Equation (|23|) may be rearranged to provide the thin 
shell's equation of motion given by 



+ V{a) = . 



The potential is given by 
V{a) = F{a) 



ms{a) 


2 


'aG{ay 


2a 




ms{a)_ 



(28) 



(29) 



where, for notational convenience, the factors F(a) and 
G{a) are defined as 

s ?7i(a) + M , „, , M — m(a) , 

F{a) ^ 1 and G(a) = ^-^ . (30) 

a a 

Linearizing around a stable solution situated at oq, we 
consider a Taylor expansion of V{a) around ao to second 
order, given by 

Via) = V{ao) + V'{ao)ia-ao) 

+ lv"{ao){a-ao)' + 0[{a-ao)'] . (31) 

Evaluated at the static solution, at a = oq, we verify 
that V{ao) = and V^'(ao) — 0. From the condition 
y (ao) = 0, one extracts the following useful equilibrium 
relationship 



2ao 



\ TTls J \ m. 



,(32) 



which will be used in determining the master equation, 
responsible for dictating the stable equilibrium configu- 
rations. 

The solution is stable if and only if V{a) has a local 
minimum at ao and ^"(ao) > is verified. Thus, from 
the latter stability condition, one may deduce the master 
equation, given by 



m 



da 



>e. 



(33) 



by using Eq. (I26p , where 770 
simplicity, is defined by 



7;(ao) and 0, for notational 



(27) with 
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(34) 



(35) 
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Now, from the master equation we find that the stable 
equilibrium regions are dictated by the following inequal- 
ities 



Vo > ^, 

Vo < ^, 
with the definition 

n = e 



if 
if 



da 



da 
da^ 
da 



>0, 



<0, 



(36) 
(37) 

(38) 



C. Stability regions 

We now determine the stability regions dictated by 
the inequalities ([36)) - (|37)) . In the specific cases that fol- 
low, the explicit form of is extremely messy, so that 
we find it more instructive to show the stability regions 
graphically. 

For the case of interest under consideration, namely, 
/3 < 1.9, we verify that da^/da\ao < 0, so that the sta- 
bility regions are dictated by inequality (|37l) . The latter 
is shown graphically in Fig. [21 for the specific case of 
/3 = 1.0. 




-0.01 2 2.2 2.4 2.6 2.8 

a/M 



3.2 3.4 



FIG. 2: Sign of {da^ /da)\ag, for 13 = 1, which corresponds to 
the absence of horizons. 



The above analysis shows that stable configurations 
of the surface layer, located sufficiently near to where 
the event horizon is expected to form, do indeed exist. 
Therefore, considering these models, one may conclude 
that the exterior geometry of a dark energy star would 
be pract ically indistinguishable from a black hole. 




stability region 



2.5 



a 



3.5 



4.5 



FIG. 3: Stability plot for 13=1, which corresponds to the 
absence of horizons; we have defined a = a/M. 



IV. CONCLUSION 

In this work, we have found exact gravastar solutions 
in the context of noncommutative geometry, and briefly 
explored their physical properties and characteristics. 
The energy density of these geometries is a smeared and 
particle-like gravitational source, where the mass is dif- 
fused throughout a region of linear dimension ^/a due to 
the intrinsic uncertainty encoded in the coordinate com- 
mutator. 

We further explored the dynamical stability of the 
transition layer of these dark energy stars to linearized 
spherically symmetric radial perturbations about static 
equilibrium solutions. It was found that large stability 
regions do exist, which are located sufficiently close to 
where the event horizon is expected to form, so that it 
would be difficult to distinguish the exterior geometry of 
the gravastars, analyzed in this work, from a black hole. 



Considering the case of /? = 1.0, the respective stability 
regions are given by the plot depicted below the curve 
in Fig. 131 Note the existence of large stability regions 
sufficiently close to the event horizon. For this case, the 
stability regions decrease for increasing a, and increases 
again as a increases. 
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